% Matlab/Octave script% Routine to write out specific beta sets for % simple curve fitting.% SWH Feb 2012close allclear all% load dataTread    % Writting out of nn_th beta for all gamma and epsilonnn=12;etemp=zeros(nbeta,neps);ptemp=zeros(ngamma,neps);wtemp=zeros(ngamma,neps);    for jj=1:neps   glow=(jj-1)*leps+1+nn*ngamma;   ghi=glow+ngamma-1;   ptemp(:,jj)=thetap(glow:ghi);   wtemp(:,jj)=thetaw(glow:ghi);   for kk=1:nbeta      blow=(kk-1)*ngamma+1+(jj-1)*leps;      etemp(kk,jj)=thetae(blow,1);   endend% note: matlab may require uppercase ASCII & TABS (not sure!)save('..//results//e_eps.txt','etemp', '-ascii', '-tabs');save('..//results//p_eps.txt','ptemp', '-ascii', '-tabs');save('..//results//w_eps.txt','wtemp', '-ascii', '-tabs');save('..//results//vgamma.txt','vgamma', '-ascii', '-tabs');save('..//results//vbeta.txt','vbeta', '-ascii', '-tabs');save('..//results//veps.txt','veps', '-ascii', '-tabs');clear all    